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Abstract. 

This review deals with properties of very young (< 10-20 Myr) stellar 
clusters which are in the nebular phase and are embedded in photoionized 
regions (classical extragalactic HII regions (RHIIs) or starburst galaxies). 
Based on the analysis of the integrated light of these clusters at the UV 
and optical wavelengths, different techniques are discussed that allow to 
estimate in consistent the stellar content and the evolutionary state of 
the ionizing clusters. Figures illustrating this contribution can be found 



1. Introduction 

RHIIs are amongst the brightest objects in galaxies. They have sizes up to 
a few 100 pc and their Ha luminosity is larger than 10^^ erg s~^ (Kennicutt 
1984). The most luminous and brightest RHII in the Local Group is 30 Doradus 
in the LMC. It has a nebular size of 400 pc and an Ha luminosity lO^'^'^ erg 
s~^, produced by more than one thousand ionizing stars. However, most of the 
radiative and mechanical energy of the region is provided by the central compact 
cluster R136. It contains more than one hundred stars within a few parcsecs, 
having a stellar density of ~ 2 stars/pc^ (Campbell et al 1992). This cluster is 
very young, ~ 3 Myr (Vacca et al 1995). Its compactness contrasts with that 
of the second brightest RHII in the Local group, NGC 604 in M 33. As 30 Dor, 
NGC 604 is photoionized by a very young cluster (3 Myr old, Gonzalez Delgado 
& Perez 2000) that contains a large population of O and B stars. However, it 
extends in a larger area with an effective radius of 60 pc (Hunter et al 1996; 
Drissen et al 1993). Several subclumps are detected, but they have still a stellar 
density is which a factor 10 lower than R136. The stellar clusters in these 
two luminous and brightest RHIIs could represent two different modes of star 
formation: compact stellar clusters (called super star clusters SSCs) and a young 
extended population that may be responsible of the diffuse UV light detected in 
more distant star forming regions. 

Starbursts show properties (nebular size, luminosity, ...) that could be con- 
sidered as scaled up versions of RHIIs where these two modes of star formation 
also co-exist (see Maiz-Apellaniz in this volume). Meurer et al (1995) have 
found that only 20% of the UV light of a sample of optically or UV selected 
starburst galaxies is provided by SSCs, and the remaining 80% could be due to 
an extended unresolved population similar to the population in NGC 604. In 
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fact, Tremonti et al (2001) have shown that this diffuse UV hght detected in the 
starburst galaxy NGC 5253 is not scattered Ught from the SSCs but it originates 
from unresolved, more extended stellar clusters. 

Most of the RHIIs and starbursts are far-away enough that the main charac- 
teristics of the stellar clusters can be estimated only through the analysis of the 
integrated light. UV light is potentially very useful to estimate the properties 
of these stellar clusters because it represents direct emission from hot massive 
stars. Therefore, UV light provides direct evidence of the location and properties 
of the most recent stellar clusters in a star forming galaxy. On the other hand, 
these massive stars arc the main source of the ionization, mechanical heating 
and chemical enrichment of the interstellar gas. The spectrum of the ionized 
gas is dominated by nebular emission. Here, I will show how to unveil the prop- 
erties of young stellar clusters embedded in photoionized regions through the 
modelling of their UV and optical light. 

2. UV-light 

The UV spectra of RHIIs and starbursts are mainly dominated by absorption 
lines (Rosa et al 1984; Kinney et al 1993). They form in the wind and pho- 
tosphere of massive stars. However, these resonance lines can also form in the 
interstellar medium of the star forming regions. Instead of strong absorption 
lines, some RHIIs show a few or weak nebular lines (see for example NGC 2363A, 
Drissen et al 2000). An inventory of the most important absorption lines are: 

• Wind lines: The strongest one are OVI A1036, NV A 1240, SilV A1400, CIV 
A1550, Hell A 1640, NIV A1720. They show a P-Cygni profile and/or are 
shifted by ~ 1000-3000 km s~^. These lines depend on stellar luminosity, 
mass loss rate, terminal velocity and age of the stars (Leitherer & Lamers 
1991). Therefore, the shape of the line profiles in the UV integrated light 
of a stellar population is related to its content in massive stars. Thus, 
these features can be used to constrain the properties (such as age, and 
the slope and upper mass limit of the IMF) of the stellar clusters 

• Photospheric lines: they are much weaker than the wind lines. SV A1502, 
FeV A1430, Silll A1417, FeV A1360-1380, OV A1371 are produced by O 
and early B stars; and Sill A1265, 1485, Silll triplet (A1295-1300), CII 
A1335 and CII A1247, 1427 also form in late type B stars. At longer 
wavelengths, Fell A2570-2615 and Mgll A2780-2825 features also form. 
These lines are also useful to constrain the age of the stellar cluster, but 
also the metallicity. In particular blends of these lines in 1360-1380 and 
1415-1435 show a strong dependence with Z (Leitherer et al 2001). 

• The most important interstellar lines in the ~1000-1800 A spectral range 
are: Sill A1260, Ol+Sill A1302, CII A1335, SilV A1400, Sill A1526, CIV 

A1550, Fell A1608, AlII A1670. Low ionization lines are very useful to 
study the kinematics of the ionized gas because the interstellar component 
usually dominates the stellar contribution (Gonzalez Delgado et al 1998a). 
They arc also useful to derive the metallicity of the gas (Pettini et al 
2000). The high ionization interstellar lines are blended with the wind 
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lines, a careful separation between both components is necessary. However, 
when the stellar cluster is young (2-8 Myr), wind lines dominate over the 
interstellar ones. 

Evolutionary synthesis models that predict the wind and photospheric lines 
have been proposed by Mass-Hesse &; Kunth (1991), Robert et al (1993), Lei- 
therer et al (1995a), Gonzalez Delgado et al (1997b), de Mello et al (2000). 
These models have been extensively applied to date stellar clusters in RHIIs 
and Starbursts observed by lUE and HST (Table 1). Most recent results are 
commented in section 2.2. 

2.1. Estimating the masses of stellar clusters 

The mass of the stellar cluster can be estimated comparing the UV luminosity 
provided by the cluster with the predictions by evolutionary synthesis models. 
However, the mass estimated represents a lower limit to the virial mass of the 
cluster because the UV light does not account for low mass stars. On the 
other hand, these estimates depend on the extinction correction applied to the 
observed UV flux. For young clusters not heavily obscured that emit a significant 
fraction of their bolometric luminosity in the UV, the extinction can be estimated 
comparing the observed with the UV intrinsic spectral distribution predicted by 
the evolutionary synthesis models, because the UV spectral slope (Fjjv ^ is 
independent of the IMF and the star formation history. It is correlated with the 
color excess E(B-V) estimated from the Balmer decrement and with Lfir/Luv 
(Meurer et al 1997). This indicates that some UV photons are absorbed by dust 
and re-emited at FIR wavelengths. This is the prediction of a dust screen model. 
However, variations on this geometry can also explain this correlation (Chariot 
& Fall 2000). In addition, the UV luminosity depends on the extinction laws 
used to do the correction. Several attenuation laws have been proposed (Rosa 
& Benvenuti 1994; Mas-Hesse & Kunth 1999; Calzetti 1997) to be applied to 
star forming regions, but all of them are coincident in showing a 2175 A bump 
weaker than in the galatic extinction law. 

Masses estimated for ionizing stellar clusters are given in Table 1. The 
masses of the clusters in RHIIs are ~ 10"^ M0, much lower than the mass es- 
timated for starbursts (--^ 10^ M©). The latter probably correspond to several 
clusters that contribute to a larger aperture. 

2.2. General Results 

The main results of these studies can be summarized as follows: UV integrated 
light is characterized by an instantaneous burst a few Myr old, populated by 
a Salpeter IMF with stars up to M^p > 50 M0. When the integrated light is 
emitted by extended areas 100 pc), the UV spectra are equally well fitted 
by continuous star formation lasting for a few Myr. These results indicate that 
clusters form with a very small age spread. In fact, this is the case for the 
starburst He2-10 (Johnson et al 2000) in which the clusters are chained along ~ 
90 pc with a mean separation < 10 pc, and they are all 4-5 Myr old. However, 
a good exception is the case of the low metallicity RHII NGC 2363. It contains 
two clusters (called A and B) within a distance of 70 pc. The UV spectra of 
the two clusters have a very different morphology. Cluster A shows only a few 
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Table 1. Age and mass estimated for the clusters responsible of the 
photoionization in RHIIs and Starbursts 
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weak interstellar lines plus CIV A1550 in emission, but the spectrum of cluster 
B shows strong wind lines (Drissen et al 2000). These differences are due to 
an age spread in the formation of the two clusters. Cluster A is still embedded 
in dust (age < 1 Myr), however, cluster B is more evolved as indicated by the 
strong wind CIV line detected, suggesting an age of 2.5 Myr. These differences 
in the age could be due to several bursts of star formation in a time scale of 
10 Myr over a spatial scale of 400 pc, probably triggered by the pass of a small 
satellite galaxy. 

There are indications that the IMF and the global star formation processes 
are the same in metal rich clusters as they are in metal poor ones. A good 
example is the metal-rich, barred starburst NGC 3049. HST observations done 
with STIS/MAMA (FUV) indicate that most of the UV light is emitted by the 
central 2 arcsec. The wind lines detected in the spectrum indicate that the 
cluster(s) in the inner 50 pc form in an instantaneous burst 3-4 Myr ago. Even 
though the metallicity of the stars is oversolar, stars more massive than 50 M© 
form in the cluster (s). However, field extended stellar population could form 
with a different IMF than the clusters. This is the case of the low metallicity 
starburst NGC 5253. Tremonti et al (2001) have obtained STIS/MAMA (FUV) 
narrow slit spectra of 8 stellar clusters plus several inter-cluster regions of diffuse 
light. They find that the UV light of all clusters is well fitted by instantaneous 
bursts with ages between 1 and 8 Myr that follow a Salpeter IMF extending up 
to 100 M0. However, the field spectrum is better fitted by continuous star for- 
mation models with either M„p=30 M© or an IMF slope steeper than Salpeter. 
However, other more sophisticated explanations are possible involving age and 
disruption effects in time scale of ~ 10 Myr. 
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2.3. Uncertainties in the models 

The basic ingredients in the evolutionary synthesis models are: evolutionary 
stellar tracks, stellar atmospheres and stellar libraries used to predict the stellar 
parameters. Thus, the results obtained depend on these inputs to the evolu- 
tionary code. Here, I comment the limitations and improvements of the mod- 
els concerning to the UV stellar libraries, evolutionary stellar tracks and the 
parametrization of the IMF. 

UV stellar libraries Most of these evolutionary synthesis models use as input 
a stellar library built with the UV spectra of stars in the solar neighborhood. 
Then, it is very suitable to describe clusters with metallicity close to solar. 
Because, wind lines depend on mass loss rate (in consequence on metallicity), 
the predictions of these models could be no adequate for low metallicty stellar 
clusters. Recently, Leitherer et al (2001) have built a UV stellar library with 
HST observations of 03 to BO stars from the LMC and SMC. They provide new 
evolutionary synthesis models at 1/4 Z0. These models predict photospheric 
lines much weaker, as expected from lower elemental abundances. However, the 
behaviour of the wind stellar lines is more complex; while, NV A1240 and SilV 
A1400 do not scale monotonically with metallicity, CIV A1550 is significantly 
affected, showing a weaker P Cygni profile. Thus, while the wind NV and 
SilV may be equally well predicted using the solar stellar library, CIV and the 
photospheric lines are overpredicted in low metallicity clusters, inducing a wrong 
estimation of the age and of the IMF parameters. 

Evolutionary stellar tracks Massive stars are fast rotators. Rotation induces 
instabilities that produce transport of angular momentum and chemical elements 
to the outer radiative envelope. When the effect of the rotation is taken into 
account, the evolution of massive stars is significantly modified. New stellar 
tracks with rotation from the Geneva group (Maeder &: Meynet 2000) predict 
with respect to the non-rotating models: a) Massive stars in the Main Sequence 
are bluer, b) The Wolf-Rayet phase is longer and starts earlier in the evolution, 
c) The ratio of the number of Wolf-Rayet to O type stars is larger, and Wolf- 
Rayet C to Wolf-Rayet N is lower. None of the evolutionary codes developed 
until now include these new tracks. 

Stochastic effects on the IMF Evolutionary synthesis models use continuous 
functions (usually power laws) to distribute the number of stars formed in a 
cluster, and thus, they do not reproduce the discontinuous nature of star forma- 
tion. This difference between the nature and the models has an important effect, 
in particular, in stellar clusters with small number of stars. Cervino et al (2000) 
have performed Montecarlo realizations for clusters of several masses (10^ M0, 
10"^ M0 and 10^ M0) to simulate the stochastic nature of the IMF. They find 
that the properties of the stellar cluster are not affected if the mass of the cluster 
is > 10^ M0. But, for lower masses, the widths of the parameter distributions 
compared with analytical values are proportional to the mass transformed into 
stars, assumptions on IMF and age. Because the mass estimated for RHIIs is 
lower than 10^ M0, stochastic effects on the IMF is an important uncertainty 
in the predictions of the properties of these stellar clusters. 
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3. The optical-IR light 

In contrast to the UV light, the optical to infrared spectrum of RHIIs is dom- 
inated by the nebular emission lines. The interstellar gas is photoionized by 
Lyman continuum photons provided by the most massive and hot stars in the 
clusters. Gas cools down via recombination and collisionally excited nebular 
lines. However, when the stellar cluster is ^ 3-5 Myr old, broad emission lines 
(Hell A4686, CIV A5800 bumps) from Wolf-Rayets are detected in the optical 
spectra. Evolutionary synthesis models that make predictions of the equivalent 
width and luminosity of these features are in Schacrer &i Vacca (1998), Star- 
burst99 (Leitherer ct al 1999), Cerviho et al (2001). 

Other optical stellar features are the high order Balmer series and some 
Hcl lines in absorption. They form in the photosphere of O, B and A stars. 
These stellar features are normally overwhelmed by the nebular contribution if 
the cluster is in the nebular phase. However, if the stars and gas do not have the 
same spatial distribution, then, US, He, H8, H9 and Hel A4471, Hel A4026, 
Hel A3819, can be detected in absorption (Gonzalez Delgado & Perez 2000). 
The strength and profile of these lines show a strong dependency with the age 
of the stellar cluster, and they are very useful diagnostics if the cluster is older 
than ~ 5 Myr old (Gonzalez Delgado et al 1999b). In addition to the Balmer 
and Hel absorption lines, the spectra of young clusters in the post-nebular phase 
(10-20 Myr old) show also the Call triplet at 8550 A and CO A2.2 //mm band. 
Beautiful examples of clusters in this phase are most nearby SSCs (called A and 
B) in the irregular galaxy NGC 1569 (Gonzalez Delgado et al 1997a; see Gilbert 
et al in this volume and the most recent works based on HST by Origlia et al 
2001 and Maoz ct al 2001). Evolutionary synthesis models treating these lines 
include Mayya (1997), Garcia- Vargas et al (1998) and Origlia et al (1999). 

In the following sections, I discuss only the results obtained from the mod- 
elling of the nebular lines. 

3.1. Modelling of the nebular emission lines 

The modelling of the nebular lines is a degenerated problem because of the 
dependence on radiation field from the ionizing stellar cluster, the chemical 
abundance of the gas, and the density structure and geometry of the region. 
Even so, it is possible to constrain the stellar content, age and IMF of the cluster 
by coupling an evolutionary synthesis and a photoionization code. Observable 
quantities to constrain the models are: 

• Collisional excited lines: Forbidden to Balmer line ratios, as [OIII] A5007/H/3; 
r] parameter, defined as ([OIII] A5007/[OII] A3727)/([SIII] A9069+9531/[SII] 
A6717+6732) sensible to the effective temperature of the cluster (Vilchez &: 
Pagel); low ionization lines like [01] A6300/H/3; and infrared fine structure 
lines such as [Nelll]/ [Nell], [Arlll]/ [Aril] or [SIV]/[SIII] 

• Recombination lines: The ratio of ionized He to ionized H, like Hel A5876/H/3, 
depends on the ratio of He to H ionizing photons, Q(He)/Q(H); thus, it 

is sensitive to the effective temperature of the stars if Teff < 40000 K 
(Stasihska 1996). 
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• Hydrogen recombination lines equivalent width, e.g. Ew(H/3), depend on 
Q(H) and also on the optical continuum luminosity emitted by the cluster. 
It is a good indicator of the age of the stellar cluster. 

Unfortunately these observational constrains also depend on other quan- 
tities like chemical abundances, differential extinction, and the contribution of 
other more evolved stellar populations. The uncertainties associated to these 
constrains and with the models are commented in section 3.2. A detailed mod- 
elling can be done follow the illustrative flow chart proposed by Garcia- Vargas 
et al (1997). A photoionization code such as CLOUDY requires in addition to 
the spectral energy distribution (SED) from the evolutionary synthesis code the 
following inputs: a) Chemical abundance of the gas estimated from the emission 
line ratios, b) Gas electron density estimated from the [SIX] A6717/[SII] A6732 
ratio, c) Assumption about the geometry of the gas. If a sphere is assumed, 
the inner and outer radius has to be specified, that can be estimated from high 
spatial resolution Ha images, d) Total ionizing photon luminosity, that it can 
be estimated from total Ha flux or from UV continuum luminosity provided by 
the stellar clusters. 

The general results from the most recent grids of models (e.g. Stasiiiska 
et al 2001; Kewley et al 2001; Moy et al 2001; Dopita et al 2000; Bresolin et 
al 1999), tailored models at optical wavelengths (e.g. Luridiana &: Peimbert 
2001 for NGC 5461; Gonzalez Delgado & Perez 2000 for NGC 604; Stasihska & 
Schaerer 1999 for IZwl8; Gonzalez Delgado et al 1999a for the nuclear starburst 
of NGC 7714; Luridiana et al 1999 for NGC 2363; Garcia- Vargas et al 1997 for 
the RHIIs of NGC 7714) and at infrared wavelengths (e.g. Colbert et al 1999 
for Arp 299; Forster-Schreiber et al 1999 for M82; Schaerer & Stasihska 1999 
for NGC 5253 and IIZw 40) can be summarized as: a) Stellar clusters are very 
young and well described by an instantaneous burst, b) Clusters form following 
a Salpetcr IMF, but in high mctalhcity RHIIs M„p is < 40-30 M0. It seems to 
contradict the results obtained at the UV light and at the optical based on the 
detection of Wolf-Rayet features in high metallicity starbursts (Schaerer et al 
2000). 

3.2. Critical points and uncertainties in the modelling 

Most of these models fail to reproduce the Ew(H/3) and some lines ratios like 
[OIII] A4363/H/3 and [01] A6300/H/3 (sec e.g. Stasihska 2000). However, the 
input to the photoionization models (e.g. the hypothesis that the nebula is ra- 
diation bounded; constant electron density along the RHII) and the ingredients 
used in the evolutionary synthesis models (in addition to those mentioned in sec- 
tion 2.3, the stellar atmospheres assumed to built the SED) have an important 
impact on the results of the photoionization models. Changes in these ingre- 
dients help to explain the divergences between the observed quantities and the 
predictions. 

Stellar atmospheres It is not clear how realistic is the hardness of the SED pre- 
dicted by evolutionary synthesis models, in particular in the Wolf-Rayet phase. 
Most of the evolutionary codes use the Kurucz (1993) stellar atmospheres. How- 
ever, Stasihska et al (2001) have produced a grid of photoionization models using 
the COSTAR stellar atmospheres (Schaerer & de Koter 1997) and they find an 
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important impact on Q(He)/Q(H), produced by the effect of the wind of O main 
sequence stars in the models. Starburst99 (Leitherer et al 1999) which is op- 
timized for predictions for young clusters, uses the non-blanketed Schmutz et 
al (1992) model atmospheres for stars with strong winds. These models pre- 
dict a large increase in the nebular excitation at the onset of the Wolf-Rayet 
phase due to a large energy above the Hell ionization limit. In consequence, 
high-metallicity models predict much higher excitation than observed in metal- 
rich RHIIs (Bresohn et al 1999: Dopita et al 2000; Kewley et al 2001) due to a 
large number of Wolf-Rayet stars. It is expected that line blanketed extended 
atmosphere models should produced a softer spectrum above 4 Ryd but it is 
not clear how the SED could be from 1 to 4 Ryd (Schaerer 2001). These line 
blanketed stellar atmospheres should predict more similar emission line ratios to 
those observed in high metallicity RHIIs. Note that the present models predict 
older ages and a lower number of massive stars, thus, a steeper IMF and a lower 
Mup, than the real values in metal-rich stellar clusters. This may be the reason 
of the discrepancy between the results obtained analyzing the nebular lines and 
the the modelling of UV wind lines in metal rich clusters. 

Density structure Most models assume that the electron density is constant. 
However, there are evidences that RHIIs can have a density structure (Perez 
et al 2001; Castaheda et al 1992). Luridiana & Peimbert (2001) have modelled 
the RHII NGC 5461 in MlOl. They find that an asymmetric nebula with a 
gaussian electron density distribution powered by a young cluster ~3 Myr old 
reproduces most of the constraints. They find that the results strongly depend 
on the assumed density law. If it is constant, then the models overestimate 
the hardness of the ionizing field, and in consequence, it gives an erroneous 
determination of the age and IMF of the stellar cluster. 

Radiation bounded hypothesis Most of the photoionization models assume that 
the nebula is radiation bounded. Consistent models using the UV light and 
nebular emission lines (e.g. Gonzalez Delgado et al 1999a) find that the mass 
of the cluster estimated from the UV continuum is in agreement with the value 
estimated with the total Balmer recombination lines (e.g. H/3). This result 
indicates that the leakage of ionizing photons is very small and, in general, that 
the radiation bounded hypothesis is correct. However, a small leakage in some 
direction is possible (Leitherer et al 1995b; Heckman et al 2001). Photon leakage 
affects preferentially the highest frequencies. These photons reach further into 
the gas, and filaments or diffuse gas can absorb them. Thus, larger values of [01] 
A6300/H/3 can be predicted by photoionization models if there is some leakage 
of ionizing photons that are absorbed by filaments (e.g. in IZwlS Stasiiiska & 
Schaerer 1999). This result rules out shocks as additional heating source that is 
able to explain larger values of [01] A6300/H/3. 

Photon leakage has also been proposed to explain the low values of Ew(H/3) 
compared with the prediction of the models. Moy et al (2001) compute pho- 
toionization models for several covering factors (cf). However, to explain the 
Ew(H/3) observed in HII galaxies a cf=0.1 is required. This means that 90% 
of the ionizing photons have to escape from the nebula. Thus, leakage of ion- 
izing photons is not a good explanation to the low values of Ew(H/3) observed 
in HII and starburst galaxies. Other alternatives, such as differential extinction 
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between gas and stars or the contribution to the contimium of an underlying old 
stellar population have been suggested (Raiman et al 2000). 

4. Conclusions 

Very general conclusions from the studies discussed here are: 

• The modelling of the wind and photospheric UV lines and the nebular 
Optical-IR lines are powerful tools to estimate the stellar content, IMF 
and age of embedded stellar clusters in RHIIs. However, we have to be 
aware that the results depend on the ingredients of the stellar evolutionary 
tracks and stellar atmosphere models, in particular on Wolf-Rayet stars. 

• Stellar clusters arc well described by an instantaneous burst of a few Myr 
old. The stellar masses distributed with a Salpeter IMF with M^p > 50 
M0. Lower M„p estimated in high metallicity RHIIs could be related to 
the use of inadequate stellar atmosphere models to compute the SED. 

Acknowledgments. I thank to the organizers, in particular to Eva Grebel 

and Dough Geisler for inviting me to give this talk and for their financial support. 
I am also very grateful to Christy Tremonti for sending me a draft of her paper on 
NGC 5253 in advance of publication, and to Grazyna Stasiiiska, Daniel Schaerer 
and Enrique Perez for their comments from a thorough reading of the paper. 

References 

Bresolin, F., Kennicut, R.C., Garnett, D.R., 1999, ApJ, 510, 104 

Calzetti D., 1997, AJ, 113, 162 
Campbell, A., et al, 1992, AJ, 104, 1721 

Castaheda, H., Vflchez, J.M., Copetti, M., 1992, A&A, 260, 370 

Cervino, M. Luridiana, V., Castander, F.J., 2000, A&A, 360, 5 

Cervino, M., Mas-Hesse, J.M., Kunth, D., 2001, A&A, submitted 

Chariot, S., Fall, M., 2000, ApJ, 539, 718 

Colbert, J.W., et al, 1999, ApJ, 511,521 

Conti, R, Leitherer, C, Vacca, W., 1997, ApJ, ApJ, 461, 87 

de Mello, D.F., Leitherer, C, Heckman, T.M., 1999, ApJ, 530, 251 

Dopita, M. et al, 2000, ApJ, 542 

Drissen L., Moffat, A.F.J. , Shara, M., 1993, AJ, 105, 1400 
Drisscn L., Roy J-R., Robert C, Dcvost D., Doyon R. 2000, AJ, 119, 688 
Forster-Schreiber, N., 1998, PhD thesis, Ludwig-Maximilian Universitat 
Garcia- Vargas M.L., Gonzalez-Delgado R.M., et al., 1997, ApJ, 478, 112 
Garcia- Vargas M.L., Molla, M., Bressan, S., 1998, A&AS, 130, 513 
Gonzalez Delgado R.M, Leitherer, C, et al, 1997a, ApJ, 483, 705 
Gonzalez Delgado R.M, Leitherer, C, Heckman, T., 1997b, ApJ, 489, 601 
Gonzalez Delgado R.M, Leitherer, C, et al, 1998a, ApJ, 495, 698 



10 



Gonzalez Delgado 



Gonzalez Delgado R.M, Heckman, T., et al, 1998b, ApJ, 505, 174 
Gonzalez Delgado R.M, Garcia- Vargas M.L., et al., 1999a, ApJ, 513, 707 
Gonzalez Delgado R.M, Leitherer C., Heckman T., 1999b, ApJS, 125, 489 
Gonzalez Delgado R.M, Perez, E., 2000, MNRAS, 317, 64 

Gonzalez Delgado R.M, Leitherer C., 2001, in "Starburst near and far", in press 
Heckman, T., et al, 2001, ApJ, |astro-ph/01050T2 

Hunter D.A., Baum W.A., O'Neil E.J., Lynds R., 1996, ApJ, 456, 174 
Johnson, K. et al, 2000, ApJ, 120, 1273 
Kennicutt, R., 1984, ApJ, 287, 116 
Kewley, L.J., et al, 2001, ApJS, 132, 37 
Kinney, A., et al, 1993, ApJS, 86, 5 

Kurucz, R. L., 1993, CD-ROM 13, ATLAS9 Stellar Atmosphere Programs and 

2 km/s Grid (Cambridge: Smithsoniam Astrophys. Obs.) 
Leitherer C, Lamers, H., 1991, ApJ, 373, 89 
Leitherer C, Robert C, Heckman T.M., 1995a, ApJS, 99, 173 
Leitherer C, et al, 1995b, ApJ, 454, L19 
Leitherer C, Vacca W.D., et al, 1996, ApJ, 465, 717 
Leitherer C, Schaerer D., et al, 1999, ApJS, 123, 3 
Leitherer C, et al, 2001, ApJ, 550, 724 

Luridiana V., Peimbert M., Leitherer C, 1999, ApJ, 527, 110 
Luridiana V., Peimbert M., 2001, ApJ, |astro-ph/0102128 



Maeder, A., Meynet, G., 2000, ARA&A, 38, 143 



Maoz, D., Ho, L.C., Sternberg, A., 2001, ApJ, |astro-ph/0105565 



Mas-Hesse J.M., Kunth D., 1991, A&AS, 88, 399 

Mas-Hesse J.M., Kunth D., 1999, A&A, 349, 765 

Mayya, Y.D., 1998, ApJ, 482, L149 

Meurer, G., et al, 1995, AJ, 

Meurer, G., et al, 1997, AJ, 110, 2665 

Moy, E. Rocca-Volmerange, B., Fioc, M., 2001, A&A, 365, 347 
Origlia, L. et al, 1999, ApJ, 514, 960 
Origlia, L. et al, 2001, ApJ, |astro-ph/0105T95 

Perez, E., Gonzalez Delgado, R.M., Vilchez, J.M., 2001, in "The Evolution of 

galaxies. I. Observational clues", in press 
Pettini, M. et al., 2000, ApJ, 528, 96 
Raiman, D. et al, 2000, MNRAS, 316, 559 
Robert C, Leitherer C, Heckman T.M., 1993, ApJ, 418, 749 
Rosa M., Joubert M., Benvenuti P., 1984, A&AS, 57, 361 
Rosa M., Benvenuti P., 1994, A&A, 291, 1 

Schmutz W., Leitherer C, Gruenwald R.B., 1992, PASP, 104, 1164 
Schaerer D, 2001, A.S.P. Conf. Serie, 221, 99 
Schaerer D, de Koter, A., 1997, A&A, 322, 598 



Clusters in RHIIs 



11 



Schaerer D, Guseva, N.G., et al, 2000, A&A, 362, 53 
Schaerer D, Stasinska G., 1999, A&A, 345, 17 
Schaerer D, Vacca, W., 1998, ApJ, 497, 618 
Stasinska G., 2000, New Astronomy Review, 44, 275 
Stasinska G., 1996; A.S.P. Conf. Ser. 98, 232 
Stasinska G., Schaerer D., 1999, A&AS351, 72 

Stasinska G., Schaerer D., Leitherer, C., 2001, A&A, |astro-ph/010222l 
Tremonti, C., et al, 2001, ApJ, |astro-ph/0103432| 



Vacca W.D., Robert C., Leitherer C., Conti P.S., 1995, ApJ, 444, 647 
Vilchez J.M., Pagel B.E.J. , 1988, MNRAS, 231, 257 
Whitmore, B. et al, 1999, ApJ, AJ, 118, 1551 



Discussion 



Sakhihov: Did you consider a case of continuous star formation? How do you 
account for the fraction of the Lyman continuum photons which do not con- 
tribute to the ionization process, missed leakage? 

Gonzalez Delgado: UV light may be equally well fitted by instantaneous or 
continuous star formation lasting for only 3-5 Myr. However, the nebular emis- 
sion lines of RHIIs are better modelled by instantaneous bursts. Continuous 
star formation models predict larger excitation and Ew(H/3) than the values 
observed. Leakage of ionizing photons can be accounted computing photoion- 
ization models with covering factor (cf) less than 1. Moderate changes (cf lower 
but close to 1) do not alter significantly the strongest emission line ratios. Lower 
values of cf are not realistic because the mass of the clusters estimated from the 
nebular recombination lines is in good agreement with estimations done using 
the UV continuum luminosity. 

D'Antona: Do the upper mass limits you gave for a few systems include 
the new rotating stellar models of the Geneva group? 

Gonzalez Delgado: Many of the results presented here were analyzed with 
Starburst99. This code has implemented the new stellar tracks of the Geneva 
group without rotation (Meynet et al 1994; Schaerer et al 1993 and Charbonnel 
et al 1993; Schaller et al 1992). To my knowledge no evolutionary synthesis code 
has yet implemented the new stellar tracks with rotation. 

Elmegreen: Your suggestion of two modes of star formation is consistent 
with optical observations, but in the infrared, most embedded stars in the solar 
neighborhood are forming in clusters, and in this sense, there is only one mode: 
it is all clustered. So what you see optically, is perhaps more a difference in 
survival and dispersal of clusters than in modes of formation of stars. 
Gonzalez Delgado: It could be true that diffuse UV emission detected in star- 
bursts could be produced by stellar clusters that have been disrupted in a time 
scale of 10-20 Myr. Tremonti et al (2001) have proposed this interpretation for 
the starburst galaxy NGC 5253. 



